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Particles embedded in a fluctuating interface experience forces and torques mediated by the 
deformations and by the thermal fluctuations of the medium. Considering a system of two cylinders 
bound to a fluid membrane we show that the entropic contribution enhances the curvature-mediated 
repulsion between the two cylinders. This is contrary to the usual attractive Casimir force in the 
absence of curvature-mediated interactions. For a large distance between the cylinders, we retrieve 
the renormalization of the surface tension of a flat membrane due to thermal fluctuations. 

PACS numbers: 87.16. dj, 87.10.Pq 

I. INTRODUCTION 

Particles bound to an interface may interact via forces and torques of two distinct physical origins. One contribution 
to these so-called interface-mediated interactions is purely geometric and results from the deformations caused by the 
particles. Since the interface is also a thermally fluctuating medium, embedded particles may also interact through 
a fluctuation-induced interaction. The associated entropic force is an example of the more general phenomenon of 
Casimir forces between objects placed in a fluctuating medium. In its original formulation, two uncharged conducting 
plates were predicted to attract each other due to the quantum electromagnetic fluctuations of the vacuum [l| . In a 
soft matter context, fluctuation-mediated forces were, for instance, studied for objects immersed in a fluid near its 
critical point or attached to a fluid interface 0-0] • 

Interface-mediated forces have also received intense attention recently due to their possible relevance in biological 
processes: membrane-mediated interactions could aid cooperation of proteins in the biological membrane and comple- 
ment the effects of direct van der Waals' or electrostatic forces Theoretical studies of this problem have typically 
considered particles on quasi-planar fluid membranes [t^-[l|| neglecting the intrinsic nonlincarity of the underlying 
(ground state) shape equation. Some systems, especially those with a symmetry, have been studied on a nonlinear 
level without taking into account any fluctuations fl7| . 

In this paper, we investigate interface-mediated interactions in their full generality on a curved geometry including 
entropic contributions for the specific problem of two parallel cylinders bound to the same side of a membrane. The 
ground state of this problem and thus the forces at zero temperature induced by the membrane were studied in 
Refs. [l||, Gil y i a stress and torque tensors and in Refs. [l9|, [2^ via energy minimization. The method employed 
here to include thermal fluctuations is based on the calculation of the free energy of the system in a semi-classical 
approximation, where Gaussian fluctuations around the curved ground state are computed. To this end we introduce a 
new paramctrization for the fluctuation variables which is possible due to the translational symmetry of the membrane. 
The force can then be obtained by deriving the free energy with respect to the distance between the cylinders. 

II. THE MODEL 

We first start by exposing the problem and shortly retrieve the ground state configuration which will be the starting 
point for the computation of the thermal fluctuations. Consider two identical cylinders of length L and radius R bound 
to one side of the membrane, parallel to the y axis and separated by a distance d (see Fig. Q]). In the limit of large 
L/R, boundary effects at the ends of the cylinders can be neglected and the profile can be decomposed into the 
following parts: an inner section between the cylinders, two outer sections that become flat for x — > ±oo, and two 
bound sections in which the cylinder and the membrane are in contact with each other. The contact area is given 
by a c RL where a c is the wrapping angle (see Fig. [IJ. The value of a c depends on the physical situation considered: 
i) the cylinder either has a finite adhesion energy w per area so that a c is determined via an adhesion balance at 
the contact lines or ii) only a given part of the cylinder surface is adhering strongly to the membrane so that a c is 
fixed. The ground state of case i), which was studied numerically in detail in Ref. [2Cj], displays a phase diagram 
which is more complicated than in case ii). In order to avoid complications already at the ground state level we will 
thus focus on case ii) in the following by setting a c = a Q + ajj = const, where a /; is the contact angle between the 
cylinder and the outer/inner membrane. The shape of the bound parts is prescribed by the geometry of the attached 
cylinder, whereas the profiles of the free membrane sections are determined by solving the nonlinear shape equation 
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FIG. 1: (Color online) Two parallel cylinders on the same side of a fluid membrane with fixed wrapping angle a c = ai + a = 
120°. The membrane fluctuates around its ground state profile Xo- The corresponding fluctuation variable u gives the distance 
between Xo and the actual profile X measured along the normal n of the ground state membrane (see inset). 



which results from the minimization of the Helfrich Hamiltonian 12 lL \2 



H=ldA(a+^), (1) 
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where £ is the surface of the free membrane and dA is the infinitesimal area element. In this functional, a denotes the 
surface tension, k the bending rigidity, and K the local curvature of the membrane. At zero temperature the profile 
obeys translational symmetry along the y axis. It is thus convenient to introduce the angle-arc length paramctrization 
ip(s) where s is the arc length and tp the angle between the x axis and the tangent to the profile. In this paramctrization 
the curvature is given by K = ±dip/ds [+ for the inner region and — for the outer ones (see again Fig. QJ]. The shape 
equation of the surface can be written as X 2 d 2 ip/ds 2 — risuiip = with A = \J n/cr the reference length scale. The 
dimensionless quantity 77 is defined as 77 = f x /c, where f x is the force per length L of the cylinder at every point of 
the membrane (which is constant and horizontal on each membrane section). Using the stress tensor (l9j |. the zero 
temperature force on the left cylinder is given by the simple expression F^l = a(r\ — 1)L. The outer section exercises 
a pulling force — aL on the left cylinder, the force of the inner section is ar/L. Later on, we will see how the values 
of the forces will be renormalized by thermal fluctuations. For the total force without fluctuations one only has to 
determine the value of r\ which is an implicit function of d. To do so, one first solves the shape equation which admits 
different solutions V'in(s) and V'out(s) (expressed in terms of elliptic Jacobi functions) corresponding to the inner and 
outer sections and depending on the boundary conditions at each cylinder. The value of 77 in the inner section for 
any given a\ is determined implicitly by the requirement tpo = V'in(so) = cc; where sq is the arc length between the 
mid-line and the contact point on the cylinder of the inner membrane. Then sq is also implicitly determined by the 
relation 



rs 

Rsinon = / dscos-0. (2) 



The torque balance equation at equilibrium, 

R 

Ky - K Q - — (ti cos ot\ - cos a Q ) = , (3) 

where K /{ is the contact curvature in the outer/inner region, fixes the individual values a Q and aj. Solving the torque 
equation, values of 77, a , a.\ for a given d can now be determined numerically. For the case under consideration 77 < 1 
is an increasing function with the distance d so that the cylinders always repel each other This justifies that 

we have restricted our discussion to parallel cylinders. Indeed, every deviation from parallelism would directly be 
compensated by a counteracting torque. 



III. THERMAL FLUCTUATIONS 



A. The fluctuation operator H 



With the knowledge of the zero temperature profile, it is now possible to compute the entropic force. We first 
focus on the inner section and set ip( s ) = V'in(s). The position vector including fluctuations can then be written as 
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X(s, y) = Xq(s) + ( — u sin ip,0,u cos ip) where u(s,y) is the membrane fluctuation in the normal direction and Xq(s) 
the position vector of the zero temperature profile (see inset of Fig. [TJ). 

For small temperature the Helfrich Hamiltonian (JTJ) can be expanded H = Hq + to second order in u, where 
Hq is the ground state energy. The contribution V(d) of the fluctuations to the free energy is then given by: 



V{d) =(3- 1 \tl / Due 



(2) 



(4) 



where j3 1 = fc^T. The fluctuation operator can be determined by expressing the Helfrich Hamiltonian ([1} with 
the help of the parametrization X(s,y). One obtains (see appendix IX| 
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where u was assumed to satisfy periodic boundary conditions. The domain of integration is — L/2 < y < L/2 and 
—so < s < so. The thermal contribution to the force of the inner section is in principal given by = dV(d)/dd. 
However, thermal fluctuations also induce a rotation of the cylinders to maintain the torque balance. For small 
membrane curvatures the actual values of rj as well as Sq and tpo differ only slightly from their zero temperature 
values. Solving the arc length and torque equations @ and ([3]) for small deviations Srj, 5sq, Sipo, one can see that the 
inner thermal force must be corrected by a prefactor, i.e., Z(ip , 7 ?)^]n- ^ turns out that Z(ip , v) does not vary much 
from unity and is thus disregarded here. For a large curvature of the inner membrane this approximation would in 
principle break down even though a general change of the behavior is not expected. 

The computation of Eq. (jj) for every value of separation d is very difficult as has no known eigenvalues and 
cigenfunctions. To circumvent this problem we focus on the two limiting cases, the quasi-flat and the highly curved 
regime and propose an interpolating formula for intermediate separations. 



B. The quasi-flat regime 



Let us first consider the quasi- flat regime, i.e., the regime of very large d/X for which the membrane can be considered 
as flat except at the cylinders. In this case 77 ~ 1 and dipo/dd w 0. In Eq. Q u can be expanded in Fourier modes 



, cxp(i7rns/so) expiiirqy / L) with n and q two integers. An implicit cut-off A of the order of the 
inverse of the membrane thickness (a ^ 5nm) is assumed. The number of modes along the y -/s-direction are given 
respectively by N = 2ttAL and M = 4ttAsq (2so ~ d being the arc length of the inner part). Since the field u(s, y) is 
dimensional, the measure of the partition function is Dm = Y[ y s /^~ 1 dit(s, y) with /1 an arbitrary length scale which 
disappears from the expression of the force. 

To compute the energy contribution Q in the quasi-flat regime, we decompose Eq. (JSJ in two parts: 
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dsdy 



and rewrite in Fourier components according to this decomposition: 

u (— m) [G 1 {— m, n, q) + X(— m, n, q)] u (n) , 



(6) 



(7) 



where the u(n) are the Fourier components of u, with u{— m) = u{m). The terms G 1 and X correspond to the 
Fourier transform of the two terms arising in the decomposition (j6|). Namely, one has: 



G ^ ( — m, n, q) = S„ 
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which is the diagonal propagator of the flat case and 



X(-m,n,q) = - 



ds 

5ir 2 ( ran q 2 
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the non-diagonal matrix due to curvature corrections. These results allow to compute the path integral 
help of the usual formula for the integral of a quadratic weight: 
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The first term of expression (|10j) is just the free energy of the flat case: 
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The second term is the perturbation correction (with G(n, q) := G(—n,n,q)): 
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X(-m, n 2 , 9)G(n 2 , g)X(-n 2 , n 3 , q)G(n 3 , q) 



...X(-nk-i,nk,q)G(nk,q)- (12) 

A careful inspection shows that to the leading order 1/d the sum in Eq. (|12[) is dominated by terms where the 
propagator G is singular, that is at (n, g) ~ 0. As the term X(— m, n, (?) is regular at the origin, we can approximate 
the series Eq. (fT^j) by keeping only the contributions of the form X(0, 0, 0)G(n, q). This dominant contribution can be 
resummed as: 



E— 

fc=i 



fe-i 



Tr (GX) fe ~ \ ^ln [1 + E X (°' °> ?)] 



2ft 



yielding a correction to the force: 
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In the limit of large d and L the sums can be approximated by integrals. A lengthy but straightforward calculation, 
which is sketched in appendix |B] leads to the following expression of the thermal force of the inner section 

0F* . 
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with A = 



f 1 ^ 

JO l+x 



d.r 



> 0, where B 



for d/X » 1 

16C(4C 2 +3C+9) 



(15) 



and C = tan 2 (tp /A). Since is negative, 



2tt 2 (AA) 2 JO l+x 2 i + _! x -L arcta n(i) " ' 3(C+1) J 

it contributes to the repulsion between the cylinders. For d — > oo, F-^/L w — A 2 /(2n 2 /3) corresponds to the entropic 
part of the intrinsic tension as found in Ref. [24[ and denoted by r = Ti n = cr — A 2 /(27r 2 /3). Obviously, the outer 
section pulls with an intrinsic tension of opposite sign r ou t = — Tin on the left cylinder. The total thermal force 
F tt = —LA 2 AX/ ((3d) in the quasi-flat regime is thus always repulsive just as the zero temperature force F^ (see 
Fig. [2]). It should be noted that even for rj close to one the membrane is not completely flat due to the nonvanishing 
contact angle ipo- Only for ipQ = is the membrane completely flat and F s equals zero. To recover the usual attractive 
Casimir behavior one has to go beyond the 1/d expansion. 
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FIG. 2: Scaled total thermal force -F a /L = -2n/3F a /(LA 2 ) = 2irAX/d on the left cylinder in the quasi-flat regime as a 
function of separation d/X for a c — 10° (short dashes), a c — 45° (dashed-dotted line), a c = 90° (long dashes), and 120° (solid 
line) . 



C. The highly curved regime 

At small separations the membrane is highly curved if the scaled ground state force r\ is close to zero This 
highly curved regime is accessible only for large wrapping angles a c : for example, choosing the radius of the cylinders 
such that R/X = 1, then d/X > 2. If the two cylinders are in contact, i.e., d/X = 2, one obtains r\ = 0.93 for a c = 45° 
which implies that the membrane is still rather flat. For a c = 120° one is already close to the highly curved regime 
since 77 = 0.11. 



1. Change of variable 

To calculate the thermal force in this regime, consider a — zero energy cost — rigid translation of amplitude unity in 
both x and z directions. This translation of the membrane as a whole can be decomposed into normal and tangential 
components which are combinations of cos?/> and sin?/'. Now, owing to the property that tangential fluctuations leave 
H invariant (25[, it is clear that individually cos ip and sinip are zero modes of H&\ As they are also zero modes of 
the operator d 2 /dip 2 + 1, we can write in the form = \ J [ip 2 (d 2 / 'dip 2 + l)u] 2 + uH(y)udy where H(y) is a 
differential operator acting on t he variable y on ly. This form of and the fact that from the ground state we have 
ds = dip/K(ip) with K (ip) = -\/2[l — ?/ cos(i/')]/A, suggests a change of variables between s ip = ip(s). In this way, 
the s dependence is eliminated in the benefit of the new variable of integration, ip, and 
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+ m)- l /X 2 - K(ip)/2]u 2 y + u 2 yy }di>dy , (16) 

where u(s, y) has been replaced by u(ip, y) [26}. The domain of integration is now — L/2 < y < L/2 and —ipo < ip < ipo- 
This new formulation is the clue to compute the Gaussian functional integral (UJ since it relies only on the angle variable 
ip. Using the Fourier transform u(ip,y) = J2 q n v. n .qexp(iirni/j/ipo) exp(iirqy / L) with n and q two integers, we can 

write i? (2) = Sm,™,,"™,*? 5 ^": where 
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with 



where Wfc(^) := exp(— iirkip/ipo). 
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2. T/ie cut-off problem 
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(18b) 
(18c) 



In the usual Fourier decomposition u(s,y) = ^2 n u n q exp{iirns / sq) exp(i7rqy/L) an implicit cut-off A is assumed 
with JV = 2irAL and M = AttAsq the number of modes along the j/-/s-direction, respectively (see Sec. IIIIBj) . A short 
inspection shows that for a finite M in the s space there is an infinite number of modes in the vj space. Actually, an 
expansion in the tp variable of a function F(ip) has to be equal to its expansion in s space: 



\ \ v— * ^ / inmb (s) \ \—\ „ ( inns s 
F (tf ( S) ) = F n exp ^_£Uj = £ F n exp ( — 



(19) 



The relation F n = J2 P A np F p connects both kind of Fourier coefficients. It can be found by expanding tp(s) in its own 
Fourier components ipk- Actually, the expansion of the exponential of a sum of exponentials is given by the Jacobi 
Anger formula which leads to 
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with 



where Jk are Bessel functions of the first kind [27| ■ Consequently, the path integral of over 



k lf> 

N modes u n . q in the s space should involve an infinite number of modes u n . q in the ip space. To clarify this point, 



consider the Gaussian weight in the s space [29j, exp(J^ 



N<n,m<N 



u^ m S nm u n ). Wc will see in the next section that it 



can be approximated by a diagonal quadratic form in the ip space, exp(^ p u- p Sp P Up). The corresponding coefficients 



S pp are obviously given by the change of variables (j2"0j) 
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-N<n,m<N 



However, a careful analysis shows an exponential decrease of the coefficients S pp for p > M which is faster the closer 
77 to 0, but relatively slow for 77 close to 1. This relies on the fact that for small 77 the coefficient A np can be shown 



to be equal to 5 np + r] P C{-l) n - p /{n - pf (1 - 5 np ) with C 



4 £n sin yjp 



For T) = 0, A Tl 



np 
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for 



p > M; the same cut-off A can thus be implemented in s and ip space. Therefore, as long as we stay close to 77 = 0, 
we consider a constant cut-off A. This condition has to be relaxed when 77 goes to one. In Sec MI Dl we will propose 
an interpolating formula between the two regimes 77 ~ and 77 rj 1. 



3. Interaction energy for the highly curved regime 



In principle we could do a perturbative expansion of the same kind as in the quasi-flat case. A careful inspection 
shows that the small expansion parameter is 77 so that the perturbative contributions due to the off-diagonal matrix 
elements S(to, n,q) are negligible for 77 rj 0. Therefore, f3V ~ 2~ 1 TrlnD with D = S(— n, n, q) given explicitly by 



7 



S = 2ijjoL[(n 2 — l) 2 ao + (2bon 2 + do)q 2 + co<7 4 ] with the notations n = irn/i[>o and q = irq/L. The various coefficients 
are functions of d. They are given by Eqs. (fT8| with k — and do = Co/A 2 — bo/2 and can be explicitly evaluated 
in terms of Jacobi elliptic functions. Interestingly for rj = 0, we have a = 2b /X 2 = 4c /A 4 = 4/(\/2A 4 ) and 



S(m,n,q) = <5„ 
obtained in |2£ 



H(n 2 1) 



(?A) a 



(qA) /4]} which, for ipo = t, corresponds to the propagator 



for membrane tubules. In the limit of large L, the sum over q can be replaced by an integration 
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(22) 



The coefficients a, &, c± are given by a = y/ao, b = Jcq and c± = [(^/aoco ± ^o)/2] 1//2 . The calculation of the entropic 
force i 7 ^ = dV(d)/dd is now straightforward but must be done with caution. First, all coefficients a, b, c as well as ipo 
are implicit functions of d. Second, the differentiation with respect to d which is a continuous variable must be done 
at constant A even though, at first sight, the number of modes M = 2irAd is proportional to d. But M is actually the 
integer part of 2nAd and is thus insensitive to an infinitesimal change of d. Moreover, Eq. (f2"2"j) is strictly speaking only 
valid in the highly curved regime where rj rj 0. To go to larger distances one has to take into account the variation of 
the number of modes (due to our change of variable s — > tp) as well as the off-diagonal elements of S(m, n,q) which 
could be computed perturbatively. 



D. Interpolation formula for all separations 

Instead of adjusting exactly the discrete number of modes with d (which is in fact impossible), we choose to 
introduce a two parameter function g(p, a) = (pso/X) v ^ a , such that the Fourier modes n in Eq. (|2"2")l must be replaced 
by n = n/g(p, a) everywhere. This ansatz takes into account the growing number of modes with d in an approximate 
but controlled manner. The parameters p and a have to be chosen such that in the highly curved regime g(p,a) ~ 1 
whereas in approaching rj sa 1, the force should correspond to Eq. (|15[) . This ansatz has thus a second advantage: like 
a variational procedure would also do, it allows to approximate the perturbative contributions which arc very hard to 
compute. Therefore, the introduction of n is a way to interpolate between the two regimes rj sa and rj sa 1. Taking 
all this into account and replacing the ^ n by an integral, the thermal force on the left cylinder reads 

2^jf n 2 ds s dg V2 m , 

l 



[(x 2 + g 2 y)U - g 2 V 'x 2 - y 2 V] In A 



2gx 

—^{■Jx 2 — y 2 U + yV)(arctani? + + arctan_B_) 
x z 

V 

— (arctanZ? + + arctanD_) (23) 
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where x = (a d) 1/2 and y = b c . We also introduced the notations A = f±£!±£\/^[]/) , B± = V2^±gV^5 and D± = 

aV ^T ^ well as U = ^{s^^yf/ 2 - {x ^ y f/ 2 ^g)} and V = ^.{s ^(x+y)^ - ( x + y)^( So g)}. 

Asking that the large d/X limit of Eq. ([23)1 is given by Eq. (fT5| we obtain a = 5 and p = (5 exp (— 27r 2 .A^) where the 
constant of integration 5 is the only free parameter determined below. 

Inserting numerical values in Eq. ([23]) one finds that the thermal force exerted on the left cylinder by the inner 
part of the membrane is negative. It thus enhances the curvature-mediated repulsion between the cylinders. As 
the outer freely fluctuating membrane exerts a constant force F^ ut /L = A 2 /(2ir 2 (3) which does not compensate F-^ 
completely, the total thermal force F R is repulsive for all separations (see Fig. [3]). The curve of — F R shows two 
characteristic trends: at short separations d/X pa 2, the force increases with d. This is due to a fast rotation of the 
cylinders for an infinitesimal change in d (see Fig. |4} which implies that the length 2sq of the inner membrane stays 
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0.8\ 




FIG. 3: Scaled total thermal force -F a /L := -2itf3F a /(LA 2 ) = -[2tt/3^/(LA 2 ) + 1/tt] on the left cylinder as a function of 
separation d/X > 2 for R/X = 1 and a c — 10° (short dashes), a c = 45° (dashed-dotted line), a c = 90° (long dashes), and 120° 
(solid line). 




(a) (b) 

FIG. 4: (a) Inner contact angle ipo and (b) derivative of the arc length dso/dd as a function of d/X for R/X — 1 and a c = 45° 
(dashed-dotted line) and 120° (solid line). The thin solid lines correspond to the respective large d limits. 

almost unchanged, i. e., dso/dd is small. The membrane is thus more under tension and thermal fluctuations are 
strongly reduced. Since dso/dd increases with d, the force grows until it reaches a maximum value at dso/dd ~ 1/2. 
For larger separations, ipo and ds^/dd stay constant and the force — F-^ should decrease in a monotonous manner 
until it tends to the constant value LA 2 /(2ir 2 [3) in the limit d/X — >• oo. Actually, this monotonous decrease of the 
force — expected on physical grounds — can be exploited to fix the fitting parameter 5: as shown in Fig. [5l 8 has to be 
set to a value smaller than 10~ 6 . 

To summarize, the total force per length F cy \/ L = (F^j + F a )/L on the left cylinder is 

F cyl /L = -a(l - V ) + F&/L + A 2 /(2tt 2 /3) . (24) 
Since this force is always negative, there is no equilibrium position beside the limit d/X — > oo, where F cy \ — > 0. 

IV. CONCLUSION 

By developing a new approach for the computation of the free energy for a system of two cylinders bound on the 
same side of a membrane, we could evaluate the corrections caused by the thermal fluctuations to the repelling zero 
temperature force. It was found that this contribution in the section between the cylinders strongly depends on the 
membrane curvature. The calculated thermal force is always repulsive. This effect differs from the attractive Casimir 
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FIG. 5: Scaled thermal force -F?/L := -2tt/3F^/(LA 2 ) on the left cylinder as a function of separation d/X > 2 for R/X = 1, 
a c = 120° and fitting parameters 5 — 10 , 10 -2 , 10 -3 , and < 10 -6 . The first three curves reach the limit 1/tt at d/X —¥ oo 
from below. All curves with 5 < 10 -6 are identical for the given resolution of the figure and give a monotous unique force. 



force which arises from the reduction of the number of internal modes with respect to the outer ones where the ground 
state is identical everywhere. This is obviously not the case here since the zero temperature shapes of the inner and 
outer sections are different: even though the number of modes in the inner section is smaller than in the outer ones, 
their fluctuations are strongly enhanced on a curved background and thus always dominate the fluctuations of the 
outer less curved region. Non-trivial membrane geometries as the one presented here are in fact ubiquitous in nature. 
The approach of this paper is sufficiently general to calculate the physical properties of other systems with highly 
curved ground states. 
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Appendix A: Derivation of the fluctuation operator (see Eq. (f5|) ) 



To rewrite the Helfrich Hamiltonian (1) in terms of the parametrization X(s,y), one has to replace the area 
element dA and the curvature K in (1) using the functions u(s,y), ^(s), and their derivatives [23| . First, one needs 
to determine the components of the (symmetric) metric tensor g a t (a, b <G {s, y}) j30j]: 



9sy 
9yy 



dX dX 

dy ds 
fdX x 
V dy 



= u y u s , and 
1 + ul , 



(Al) 



where we have introduced the notations tp = difj/ds, u s = du/ds and u y = du/dy. This allows to calculate the area 



element dA = ^Jgdsdy with ^fg = det (g a b) = y {uip — 1) 2 (1 + u 2 y ) + u 2 s . The components of the extrinsic curvature 
tensor K a b are: 



Key = 



dX 


dn 


ds 


ds 


dX 


dn 


dy 


~ds~ 


dX 


dn 


dy 


dy 



1 



ijj + 2wt/j — u ss — u ijj — uu s ip + (uu ss — 2u s ) ip 



(-1 + Uljj) U S y - tpUyUf, 



and 



(A2) 



10 



The contraction of K a b with the metric yields the curvature K = K a tg ab where g ab is the inverse of the metric, 
i.e., g ss = g yy /g, g sy = —gsy/g, and g yy = g ss /g. Inserting the expressions for dA and K into H one identifies 

H (2) = + ^ UsU + + 3^)4 + ^2^ sU 

+ ^4^ + w ss u ra + - f )^]d S dy . (A3) 



Appendix B: Derivation of Eq. (115|) 



Let A(0) = X(0, 0, 0) and start with expression ([HI) 

1 



gJf , x V &S„*(0)G(n, g ) 
2 , 1 + 



Note that, according to their definition, ^A(O) = — ^A(O) and ^X^nG(n,q) = jEi 

continuous < 
separations: 



(Bl) 



Thus, a 



continuous approximation for the sum (which is valid for large d, given the high number of modes) yields for large 



0_ 

dd 



X)X(0)G(n,g) 



X(0) 



1 



2n 2 



(n 2 + q 2 ) (rf + cpy 



dn 



A A- 



( A2 + ? 2 ) (t) 



and by the same token: 



1 + V A(0)G(n, g ) = 1 + — X(0)d / A * _ dn ~ 1 + — arctan | | 

n d \ d A / 



As a consequence Eq. (|Blj) becomes in the continuum approximation: 



p d 2tt 



x(0)d A 2 A 



(A 2 +9 2 ) 



^ 1 + ^^ arctan 

7T q 



dg. 



A careful inspection of this expression can be performed by dividing the integration interval into three parts, 

^A 3 )' ,A 



$A,(£A») ! 



and 



It turns out that the contributions of the two first intervals 



are negligible with respect to the last one. Moreover, checking that in the range of integration considered, the 

„2 



numerator q 2 — ^A can be approximated by q 2 , we can thus write f35F for large d as: 



A X(0)d x 2 A 1 



o 1 + 



A 2 X(0)d 
7r q 



arctan ( — 



dq. 



(B2) 



Ultimatly, the evaluation of f35F requires the computation of 



A(0)=X(0,0,0) = i ' 



ds 



(B3) 



This computation relies on the saddle point solution tJj(s). Actually, for large d, 77 ~ 1 one has 



V-' 



1 — cos (ip) 



X 2 



sin 



11 



whose solution is: 



ip(s) = 2 arccos 



1 - Cexp(2^fi] 



l + Ccxp 2^ 



with C = tan 2 ($f\ . As a consequence, replacing ip(s) in Eq. ([B3]) one finds directly that X(0) = t C ^dx^c^iy^ 



Inserting this value in the expression of f35F in Eq. (|B2|) leads to the result claimed in the text 



[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948). 

[2] M. E. Fisher and P. G. de Gennes, C. R. Seances Acad. Sci., Ser. B 287, 207 (1978). 
[3] C. Hertlein et al, Nature 451, 172 (2008). 

[4] M. Kresch, The Casimir Effect in Critical Systems (World Scientific, Singapore) 1994. 

[5] H. Li and M. Kardar, Phys. Rev. Lett. 67, 3275 (1991). 

[6] H. Lehle, M. Oettel, and S. Dietrich, Europhys. Lett. 75, 174 (2006). 

[7] E. Noruzifar and M. Oettel, Phys. Rev. E 79, 051401 (2009). 

[8] B. J. Reynwar et al., Nature 447, 461 (2007). 

[9] M. Goulian, R. Bruinsma, and P. Pincus, Europhys. Lett. 22, 145 (1993); ibid. 23, 155(E) (1993). 
[10] R. Golestanian, M. Goulian, and M. Kardar, Europhys. Lett. 33, 241 (1996); Phys. Rev. E 54, 6725 (1996). 
[11] J.-M. Park and T. C. Lubensky, J. Phys. I France 6, 1217 (1996). 
[12] K. S. Kim, J. C. Neu, and G. F. Oster, Europhys. Lett. 48, 99 (1997). 
[13] P. G. Dommersnes and J.-B. Fournier, Europhys. Lett. 46, 256 (1999). 

[14] T. R. Weikl, Eur. Phys. J. E 12, 265 (2003). 

[15] C. Yolcu, I. Z. Rothstein, and M. Deserno. larXiv:1007.4760Vl. 

[16] H.-K. Lin, R. Zandi, U. Mohideen, and L. P. Prvadko. larXiv:1009.07"24V l. 

[17] P. Biscari, F. Bisi, and R. Rosso, J. Math. Biol. 45, 37 (2002); P. Biscari and F. Bisi, Eur. Phys. J. E 7, 381 (2002). P. 

Biscari and G. Napoli, Biomechan. Model Mechanobiol. 6, 297 (2007). 
[18] M. M. Miiller, M. Deserno, and J. Guven, Europhys. Lett. 69, 482 (2005); Phys. Rev. E 72, 061407 (2005). 
[19] M. M. Miiller, M. Deserno, and J. Guven, Phys. Rev. E 76, 011921 (2007). 
[20] S. Mkrtchyan, C. Ing, and J. Z. Y. Chen, Phys. Rev. E 81, 011904 (2010). 
[21] P. B. Canham, J. Theoret. Biol. 26, 61 (1970). 
[22] W. Helfrich, Z. Naturforsch. 28c, 693 (1973). 

[23] For a detailed introduction to the differential geometry of two-dimensional surfaces see: M. Do Carmo, Differential Geom- 
etry of Curves and Surfaces, (Prentice Hall, 1976); E. Kreyszig, Differential Geometry, (Dover, New York 1991). 
[24] J.-B. Fournier and C. Barbetta, Phys. Rev. Lett. 100, 078103 (2008). 
[25] R. Capovilla, J. Guven, and J. A. Santiago, J. Phys. A: Math, and Gen. 36, 6281 (2003). 
[26] The Jacobian is an infinite constant that does not affect the physical quantities and can thus be disregarded. 
[27] Handbook of Mathematical Functions, 9th ed., ed. by M. Abramowitz and I. A. Stegun (Dover, New York, 1970). 
[28] J.-B. Fournier and P. Galatola, Phys. Rev. Lett. 98, 018103 (2007). 

[29] For the sake of simplicity, we suppress the variable q which does not play any role in the following argument and write 

S nm for S(n, m, q). 
[30] Note that dX /ds = (cos^, 0, sin ip). 



